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Abstract 

Orbital magnetism in an integrable model of a multichannel ring with 



^ . long-ranged electron-electron interactions is investigated. In a noninter- 

fj ' acting multichannel system, the response to an external magnetic flux is 

the sum of many diamagnetic and paramagnetic contributions, but we 
find that for sufficiently strong correlations, the contributions of all chan- 
nels add constructively, leading to a parity (diamagnetic or paramagnetic) 
which depends only on the total number of electrons. Numerical results 
^— N ' confirm that this parity-locking efi^ect is robust with respect to subband 

(^ ' mixing due to disorder. 

^^ , The free energy F{(j)) of a metallic ring threaded by an Aharonov-Bohm flux 

r~^ ■ {hc/e)(j) is a periodic function of (/>, with period 27r [|l|. The system is said to 

^\ , be diamagnetic if F is minimal for = 0, and paramagnetic if F is minimal for 

(j) = TT (1/2 flux quantum). A purely one-dimensional (ID) ring with N spinless 
2 ! electrons is diamagnetic if N is even and paramagnetic if N is odd, independent 

P^ ■ of disorder and interactions [^ |j. This connection between the parity of N 

and the sign of the magnetic response of the system is sometimes referred to 
2 ' as Leggett's theorem [Q. The persistent current / ~ ~~{e/K)dF/d4> is a periodic 

Q \ function of (p with amplitude Iq = evp/L in a clean ID ring y|, where vp is 

O ' the Fermi velocity and L the circumference of the ring. In a ring with many 

independent channels, / is the sum of many such diamagnetic and paramagnetic 
contributions, and thus has a random sign and small amplitude [y. Leggett's 
theorem is thus generally violated in multichannel systems. 



j^ ■ In an interacting system, the parities of different channels are no longer in- 

dependent. In this article, we investigate the orbital magnetism of ID SU{M) 
fermions interacting via the potential V{x) — g/x^. We find that for 5 > 0, 
Leggett's theorem is restored for an arbitrary number of channels M . In ad- 
dition, provided the ring is sufficiently thin {kpL > 27rAf) we find that the 



magnetic response of all channels add constructively, leading to a large enhance- 
ment of the persistent current. A disordered two-channel ring with interchannel 
interactions is also investigated numerically, and shows, importantly, that the 
parity-locking effect persists even when the subbands are mixed strongly by dis- 
order. Interchannel correlations of this sort may be important to explain the 
anomalously large observed value of the persistent current in normal metal rings 

Interacting spinless electrons in a non-disordered ring with M transverse 
channels, threaded by an Aharanov-Bohm flux {hc/e)(f), may be represented by 
ID SU{M) fermions. The transverse degrees of freedom may be represented 
by an SU{M) spin variable a = 1, . . . , M. In the absence of disorder, and for 
interactions which depend only on the electrons' coordinates along the ring {i.e., 
for thin rings), the number of electrons K^r in each channel is conserved. The 
Hamiltonian of the system is 

i—l * i<j 17— 1 

where N = ^^ K^ is the total number of electrons and Eo- is the energy mini- 
mum of subband a. Units with h = m = 1 are used. The Aharanov-Bohm flux 
leads to the twisted boundary condition ||l| 

^{xiai,-- ■,{xi + L)ai,- ■■ ,xn(tn) = e*'^*(a;icri, • • ■ ,Xi(7i,-- ■ ,xn<Jn)- (2) 

For simplicity, let us consider equally spaced subbands Ecr+i — Eo- = A = Ep/M. 
The subband splitting A plays the role of an SU{M) magnetic field. As we 
shall see, the effect of repulsive interactions is to renormalize this effective field, 
causing a condensation of electrons into the lowest subband. At T = 0, the 
equilibrium persistent current is given by /(0) = —{e/h)dEo/d(f), where £^o(0) 
is the ground state energy of Eq. (||), subject to the boundary condition (||). 

The first harmonic of the persistent current may be characterized by the 
value of / at = 7r/2 (1/4 flux quantum), assuming the higher odd harmonics 
are small. For V{x) — 0, one finds 

'i<^«m-^T.i-»''"<'- (3) 

cr=l 

For M > 1, this leads to the well-known result |/| ~ AI^/'^Iq due to the 
random parities of the different channels. The system may be either diamagnetic 
or paramagnetic, depending on the channel occupancies Ka-. 

Let us next consider a model with long-range interactions: V{x) = g/d{x)'^, 
where d{x) = {L/Tr)\ sin(7rx/L)| is the chord length along the ring. This model 
was introduced and solved by Sutherland M for the case M = 1 and = 0. 
First, we show that the SU{M) model with twisted boundary conditions is 



completely integrable. Let us introduce the following generalized momentum 
operators 

7T, ^p^-iX-Y^ cot ^ \ ' Pjk , (4) 

where the operator Pjk permutes the isospins and the positions of the two 
particles j and k simultaneously, and A = y/g + 1/4 + 1/2. These generalized 
momentum operators are hermitian, tt^ = ttJ . The Hamiltonian can be expressed 
as a quadratic form in terms of the tt^, 

N M 

i? = ^ -Trf + ^ e„K„ + const. (5) 

i=l cr=l 

Introducing a set of new operators ttj = tt^ + A^ ■/ /^n Py, one can show that 
the following infinite number of physical operators commute with each other 

[InJm]=0, (6) 

where In = '^i^ii^^i)"' i with n — 0, 1, 2, . . . , oo. Furthermore, one can show 
that these operators also commute with the Hamiltonian 

[H, in] - 0. (7) 

This therefore provides a proof of the complete integrability of the system. 
Our proof of integrability follows that given by Polychronakos for the Calogero- 
Sutherland model |^. Note, however, that the operator Pij permutes the posi- 
tions and the isospins of two particles i and j, while the operator Mij of Ref. |^ 
only exchanges the positions of the two particles. 

Consider the following Jastrow wavefunction for given channel occupancies 

Wo(xiCri,X2Cr2,---,2^7VCrAr) = J_ J_ | Sm( ^-^7r)| Sm( 'T^^'^)^ \^) 

l<i<j<N 

where the isospins of the particles are symmetric under exchange. This wave- 
function satisfies periodic boundary conditions {4> = 0) if iV is odd and satisfies 
antiperiodic boundary conditions ((/> = tt) if iV is even. It also satisfies Fermi- 
Dirac statistics. Furthermore, this wavefunction is annihilated by the operators 
TTi, with i = 1,2,- ■ ■ , N: 

Tri'^o{xiai,X2(72,- ■■ ,XNcrN) ^0, (9) 

which indicates that it is the ground state of H in the subspace with fixed Ka^ . 
The absolute ground state of the system thus occurs for (p = Q when N is odd 
and for (f> — tt when N is even, regardless of the subband spacing Sa-- The 



Leggett theorem thus holds for this mteracting muhichannel system, i.e., the 
ground state is diamagnetic if N is odd and paramagnetic if N is even. 

Now, let us consider general values of the external magnetic flux. For g > 0, 
the ground state is highly degenerate in the limit L ^ oo in the absence of 
SU{M) symmetry breaking (A — 0) due to the strong repulsion of the potential 
at the origin, which prohibits particle exchange. In a finite ring, one expects 
many SU{M) level crossings |10| as a function of 4>- These states differ in energy 
by at most nfivp/L due to boundary effects; all electrons will thus be condensed 
into the lowest subband for A > nfivp/L, i.e., for kpL > 27rM, which is satisfied 
provided the ring is sufficiently thin. The ground state of the system in this 
"ferromagnetic" state has the Jastrow product form 



*({a;}) = cxp 




sm 



-1 



sm 



(10) 

for < I?!) < TT, where a = if TV is odd and a = tt if A^ is even. One readily 
verifies that \1/ is an eigenstate of Eq. (y) , has the correct symmetry, and obeys 
the twisted boundary condition (g). This state coincides with the absolute 
ground state ^o for (p = (tt) when N is odd (even) . The ground state energy 
is found to be 



6L2 2 V L , 

The corresponding persistent current is 

/(</-o/4) = (-1)^^ - {-ifMIo. (12) 

The condensation of all electrons into the lowest subband caused by the strong 
repulsive interactions thus leads to an enhancement of the typical persistent 
current in the ballistic regime by a factor of M^''^ compared to that for nonin- 
teracting electrons, given in Eq. (||). 

A peculiarity of the integrable model considered above is that the number of 
electrons in each channel is a constant of the motion. Both disorder and more 
realistic interactions which depend on the transverse coordinate will break this 
symmetry, and it is therefore important to verify that the parity-locking effect 
is not destroyed. To this end, we have considered a disordered two-channel ring, 
modeled in the tight-binding approximation, with a nearest-neighbor interchain 
interaction V included to induce interchannel correlations. The Hamiltonian is 



^ = E 



i=l 



2 . 

V (e"''/^4^C^+ia + H.C. + e^aP^c^ + TT (4^*2 + H.C.) + Vp.aPi: 



a=l 



2 , 

(13) 



where cj^, creates a spinless electron at site i of chain a, pia = c\^Cia, and Sia is a 
random number in the interval [— e/2,e/2]. The interchain hopping determines 
the subband splitting A. Fig. 1 shows the persistent current for rings with 5 
spinless electrons on 18 sites as a function of V , calculated using the Lanczos 
technique. Both the ensemble average of / for 500 systems (squares) and the 
persistent current of a typical system (solid curve) are indicated. The error bars 
indicate the standard deviation 51 = ((/^) — (/)^)^/^ over the ensemble (not the 
statistical uncertainty in the mean value). The subband splitting A — 0.8 is 
chosen so that in the absence of disorder and interactions, Ki = 3 and K2 = 2, 
leading to a large cancellation of the persistent current due to the different 
parities of the two channels. The on-site disorder e = 2 > A, Ep mixes the two 
channels, but does not lead to strong backscattering. For V^ = 0, the sign of 
I is essentially random, and (/) is slightly diamagnetic. As V is increased, (J) 
oscillates in sign and increases in magnitude, becoming strongly diamagnetic 
for large V. \{I)\ is increased by a factor of 10 as V is increased from to 20, 
while \/{P) is increased by a factor of 3. The system exhibits parity locking 
for large V. 

While the subband occupancies are no longer constants of the motion in Eq. 
(O), there is a corresponding topological invariant in the disordered system, 
namely, the number of transverse nodes in the many-body wavefunction H 
{i.e., nodes which encircle the AB flux (p). The lowest subband has no such 
nodes, while each electron in the second subband contributes one transverse 
node. In order for two electrons to pass each other as they circle the ring, such 
a transverse node must be present. As V increases, it becomes energetically 
unfavorable for electrons to approach each other, so transverse nodes in the 
many-body wavefunction will tend to be suppressed. In the strongly-correlated 
limit, all such nodes will be eliminated, leading to a state whose parity depends 
only on the total number of electrons. In such a state, the persistent currents 
of all channels add constructively, leading to a large persistent current (see 
Fig. 1). It should be emphasized that while the enhancement of the persistent 
current shown in Fig. 1 is relatively modest, a much larger enhancement would 
be expected in a system with M 3> 1, based on the above arguments. 

In conclusion, it has been shown that for sufficiently strong electronic cor- 
relations, the sign of the orbital magnetic response in a multichannel ring is 
determined solely by the parity of the total number of electrons in the system. 
The magnetic response of all channels add constructively, leading to a large 
enhancement of the persistent current. This result, which we refer to as the 
parity- locking effect, was demonstrated in an integrable model with long-range 
interactions as well as in a disordered two-channel system. It should be em- 
phasized that the parity-locking effect holds in both ballistic and disordered 
systems; it is therefore complementary to mechanisms previously proposed for 
the enhancement of the persistent current [O, O, O, n4, na , which rely on the 
competition between disorder and interactions. It is likely that both such mech- 
anisms are important to explain the anomalously large observed value p, m of 



the persistent current in disordered metallic rings. 
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Figure 1: Persistent current / = —{e/h)dEo/d(j)\ff,^^g/4 of a disordered two- 
channel ring with 5 spinless electrons on 18 sites as a function of the interchain 
interaction V. The current is given in units of /q = evp/L. The amplitude of 
the on-site disorder is £ = 2 and the subband splitting is A = 0.8. Solid curve: 
persistent current for one realization of disorder. Squares: ensemble average (/) 
for 500 systems. The error bars indicate the width SI = {{P) — {I^Y^ of the 
current distribution. Note that the persistent current is diamagnetic for large 
V, as expected for a system with N odd due to parity-locking; this is true for 
all realizations of disorder. 



